On Iterative Algorithms for Quantitative Photoacoustic Tomography in the Radiative Transport Regime

In this paper, we describe the numerical reconstruction method for quantitative photoacoustic tomography (QPAT) based on the radiative transfer equation (RTE), which models light propagation more accurately than diffusion approximation (DA). We investigate the reconstruction of absorption coefficient and/or scattering coefficient of biological tissues. Given the scattering coefficient, an improved fixed-point iterative method is proposed to retrieve the absorption coefficient for its cheap computational cost. And we prove the convergence. To retrieve two coefficients simultaneously, Barzilai-Borwein (BB) method is applied. Since the reconstruction of optical coefficients involves the solution of original and adjoint RTEs in the framework of optimization, an efficient solver with high accuracy is improved from~\cite{Gao}. Simulation experiments illustrate that the improved fixed-point iterative method and the BB method are the comparative methods for QPAT in two cases.Comment: 21 pages, 44 figure

On Relaxed Averaged Alternating Reflections (RAAR) Algorithm for Phase Retrieval from Structured Illuminations

In this paper, as opposed to the random phase masks, the structured illuminations with a pixel-dependent deterministic phase shift are considered to derandomize the model setup. The RAAR algorithm is modified to adapt to two or more diffraction patterns, and the modified RAAR algorithm operates in Fourier domain rather than space domain. The local convergence of the RAAR algorithm is proved by some eigenvalue analysis. Numerical simulations is presented to demonstrate the effectiveness and stability of the algorithm compared to the HIO (Hybrid Input-Output) method. The numerical performances show the global convergence of the RAAR in our tests.Comment: 17 pages, 26 figures, submitting to Inverse Problem